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We suggest a "minimal model" for the 3D turbulent energy spectra in super- 
fluids, based on their two-fluid description. We start from the Navier-Stokes 
equation for the normal fluid and from the coarse-grained hydrodynamic 
equation for the superfluid component ( obtained from the Euler equation for 
the superfluid velocity after averaging over the vortex lines) and introduce 
a mutual friction coupling term, proportional to the counterflow velocity, 
the average superfluid vorticity and to the temperature dependent parame- 
ter q = «/(! + a'), where a and a' denote the dimensionless parameters 
characterizing the mutual friction between quantized vortices and the normal 
component of the liquid. We then derive the energy balance equations, taking 
into account the cross-velocity correlations. We obtain all asymptotical solu- 
tions for normal and superfluid energy spectra for limiting cases of small/big 
normal to superfluid density ratio and coupling. We discuss the applicability 
of our model to superfluid He II and to ^He-B. 

PACS numbers: 67.40.Vs, 67.57.De, 47.27.Ak 



INTRODUCTION 

Quantum turbulenc^l^ containing a tangle of singly quantized vortex 
line - such as turbulence in superfluid ^He (He-II) and in the superfluid 
B-phase of ^He (^He-B) - besides of being for half a centurjEl a playground 
for low temperature physicists increasingly attracts attention of the fluid 
dynamics community. In this paper, we address the important question 
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of the turbulent energy distribution between scales (3D energy spectra) of 
such turbulence. In view of very sparse experimental data - besides indi- 
rect indications of its form in He 11^ ^ and ^He-!^ deduced from various 
decay measurements or computer simulation j^ l ^ l ^ l we recall the only direct 
measurement of the velocity spectrum based on pressure fluctuations in He 
n by Maurer and Tabelin^ - there is a clear call to tackle this important 
issue theoretically. We consider the simplest case of classically generated 
turbulence in quantum liquids such as He H or ^He-B that can be thought 
of as isothermal, homogeneous and isotropic. We work within a framework 
of the two-fluid model, building on ideas first in trodu ced by VolovilP and 
VineiP, by further developing our previous wor k I I. We stress that our 
approach does not directly apply to quantum turbulence generated in su- 
perfluid helium by the thermal counterfloW^, which is anisotropic, being 
generated thermally, by the temperature gradient in the channel. 

Above the pressure dependent transition temperature (Tx ~ 2.17 K for 
'^He and Tc ~ 1 mK for '^He) both ^He and ^He are ordinary viscous fluids 
that can be described by the Navier-Stokes equations and their turbulent 
flow is fully classical. From hydrodynamical viwepoint normal ^He and ^He 
liquids differ from each other mainly because of their very different values 
of kinematic viscosity. While liquid ^He above Tx possesses the lowest kine- 
matic viscosity, u, of all known fluids"*^^, of order 1/4 ~ 2 x 10~^ cm^/s (fifty 
times smaller then that of water at room temperature), liquid ^He at mil- 
likelvin temperature is a Fermi liquid^^ (1^3 oc T~^) with kinematic viscosity 
exceeding that of air (which is about 0.15 cm^/s) comparable with olive 
oil, of order ~ 1 cm^/s . In principle, both these liquids may become 
turbulent. 

The plan of the paper is as follows. After this short introductory Section 
we introduce the continuous minimal model for two-fluid turbulence with 
mutual friction in Section 1 and we use it to derive the turbulent energy 
spectra in Section 2. We discuss the applicability of this continuous model 
and specify the crucial role of Kelvin waves in Section 3. We conclude and 
outline the future work in Section 4. 



1. MINIMAL MODEL FOR TWO-FLUID TURBULENCE 
WITH MUTUAL FRICTION 

lA. Basic equations for the two-fluid model 

In this paper we adopt the simplest form of the two fluid model for 
superfluid ^He and ^He (see e.g. Eqs. (2.2) and (2.3) in the Donnelly's 
textboolPI) which neglects both bulk viscosity and thermal conductivity. 
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These are the Euler Eq. for the superfluid velocity Ug (with zero viscosity 
= 0) and the Navier-Stokes Eq. for the normal component (with the 
kinematic viscosity i^n = z^): 



Ps 
Pn 



' + Ms • V It, 



dt 

[dUn 

dt 



Vps = -i^ns , (la) 

+ (tin • V)U„ - Vpn - PnV/^Ui + F^s ■ (lb) 



Here p^, Ps are the densities of the normal and superfluid components, p^, 
Ps are corresponding pressures: 

Pn = —[p + Ps\us-Ua\'^], (Ic) 
P 

Ps = —[p-pn\Us-U^\'^], (Id) 
P 

Here p = Ps + Pn and F^s describes the mutual friction: 

Fns = -Ps{a'{us -u^) X u + a Us X X {us - Un)]} • (le) 

o^s = V X Its is the superfluid vorticity; a>s = u^s/^s', ol and a are dimension- 
less parameters describing the mutual friction between superfluid and normal 
components of the liquid mediated by quantized vortices which transfer mo- 
menta from the superfluid to the normal subsystem and vice versa. For the 
flow with vortices locally aligned with each other these parameters enter 
the reactive and dissipative forces acting on a vortex line as it moves with 
respect to the normal component. Here we consider o! and a as phenomeno- 
logical parameters, assuming the general case where quantized vortices are 
not aligned locally and thus the bare parameters are renormalized. 

Following Eef llll we approximate the "mutual friction term" i^ns, Eq. (|lej) . 
as follows: 

Fns = Q0o{Us - Wn) ■ (H) 

Here Q ~ aps and 0o is some characteristic superfluid vorticity 



00 = VW> • (Ig) 

Approximation ((TH) accounts for the fact that the vorticity in developed 
turbulence usually is dominated by the smallest eddies in the system of scale 
77 with the largest characteristic wave- vector A;^ ~ 1/rj. These eddies have the 
smallest turnover time Tmin that is of the order of their decorrelation time. 
On the contrary, the main contribution to the velocity in the equation for 
the dissipation of the fc-eddies with intermediate WaVG- vectors k ^maxi 
is dominated by the A:'-eddies with k' ^ k. Because the turnover time of 
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these eddies t^i » Tmin, approximation (QJ) adopts self- averaging of Eq. ((Te|) 
on time intervals of interest (rmin t <^ Ty) and thus vorticity can be 
considered as almost uncorrelated with the velocity u which is a dynamical 
variable. Actually, Eq. ((T^ is the mean field approximation that neglects 
the fluctuating part of vorticity and replaces it by by its mean value. 



IB. Model for the cross- correlation {u^ ■ u^) 

Our first goal in this paper is to formulate the energy balance equations 
for the one-dimensional energy spectra Ea^k) for the normal and the su- 
perfluid subsystems, related with the simultaneous, same-point correlations 
{ui ■ Uj) as follows: 

j dkEijik) = Eij = 2 {ui ■ uj) . (2) 

a. Some definitions and known relationships. To find the cross- 
correlation {ui ■ Uj) we need to recall some definitions and relationships, 
well known in statistical physics. The first one is the Fourier transform: 

Ui{k,0,t) = J ^^Ui(r,t)exp[-i(fc-r-0t)], (3) 

and the definition of the two-point, different-time cross-correlation functions 
Fij{k,0): 

{ui{k,0)-u*{k',0')) (4a) 
= {2T:f 5{k - k')5{0 - 0')Fij{k,0) . 

Next we need to define simultaneous two-point (cross) correlators, Fij{k), in 
fe-r epr esent at ion : 

{ui{k,t) ■ u*{k',t)) = {2Tif5{k - k')Fij{k) . (4b) 

The frequency integral of Fij{k,0) produces Fij{k): 

2-nFij{k) = j d0Fijik,0) . (4c) 

Further fc-integration gives same-point correlations, which has a sense of 
(twice) kinetic energy density per unite mass of the normal (for i = j = n) 
or the superfluid (for i = j = s) kinetic energy density per unite mass: 

(27r)-3 j SkF,.j{k) = {ui ■Uj)= 2Eij . (4d) 
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Our approach is formulated in terms of one-dimensional density of Eij 
in fc-space, 

Eij{k) = eFij{k)/{27r)\ (4e) 

defined such that 

'' dkEij{k) = Eij . (4f) 



Finally, we define a scalar version of the velocity Green's (response) function: 
5ui{k, 0) 



<5/,(fc',0') 



(27r)^,5(fe - k') (5(0 - 0') Gij{k, 0) . (5) 



To determine (tij • Uj) we adopt a few more or less justified assumptions, 
discussed below. 



b. Fluctuation-dissipation approximation for the diagonal tur- 
bulent correlation function. In the thermodynamical equilibrium, the 

different-time (cross) correlation functions (in the 0-representation) are re- 
lated with the Greens' function by the mean of the fluctuation-dissipation 
theorem, which in hydrodynamics reads: 

Fu{k,0) = 2Tlui[Gii{k,0)], (6) 

where T is the temperature of the system. Having in mind that the ther- 
modynamical equilibrium is the equipartition of the energy with T/2, being 
the energy per degree of freedom (which in our notation is Fii{k)/2), we can 
write in equilibrium: 

Fii{k,0) = 2Fu{k)lm[Gu{k,0)] . (7) 

Our conjecture is that this equation is approximately (on a semi-qualitative 
level) valid also in the flux-equilibrium state, i.e., in the fully developed 
turbulence. 

c. One-pole approximation for the diagonal Greens functions. 
Generally speaking, the Greens function in the developed turbulence is very 
involved function of k and 0. The only facts that we know for sure are: 

1. Gii{k,0) is analytical in the lower half-plane. 

/oo 
Gii{k, 0) d0 = TT , lim Gii{k, 0) = I/0 
-00 ^-"^ 
3. At given k, Gii{k,0) has a characteristic width (in 0), of about the 

total damping frequency in the system, i.e., Vik^ + 7i(fc) -|- qi0o- 
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The simplest analytical from, which satisfies these requirements is the 
so-called one-pole approximation 

^) ^ -r ,2 / ^ 1 ' (8a) 

- I [Uik^ +^i{K) + qi0o\ 

widely used in theoretical physics in general and in the theory of turbulence 
in particular. In Eq. Q Jii^k) is the eddy-decorrelation frequency, which 
is the same as the eddy-turnover frequency. The Kolmogorov 1941 (K41) 
approximation for this object reads: 

j^{k)^Ciei{k)^/''k^/\ (8b) 

where ei{k) is the energy flux in the i-subsystem. The well known effective 
turbulent viscosity v-Yik) is related with ^{k) in the following way: z^t(^) = 

-f{k)/e oc k-^/^. 

d. Approximation of the Gaussian statistics. In the limit pn 3> Ps 
the turbulent velocity w„ (or Us in the opposite limiting case, ps ^ Pn) 
can be considered as a given one, independent of Ug (or for ps ^ p„). 
Having in mind that in practice we are dealing with moderate extend of the 
inertial interval, and that the statistics of turbulence in the energy contained 
interval is very close to the Gaussian one, we can approximate the statistics 
of turbulent field w„ (or Ug, for ps » Pn) as Gaussian. 

e. Here a frog jumps. As we suggested, in the limit pn 3> Ps we can 
consider the cross- velocity term qs0oUn = fs in the RHS of Eq. ^ for Us as 
a Gaussian random force fs and compute the cross-correlation (for simplicity 
in the scalar version), (usfs), using so called Gaussian integration by parts: 
{usf*) = {6us/Sfs) (/s/s*)- With /s = qs0oun this gives for the velocity cross- 
correlation: 

Fsn{k,0) = qs0oFnn{k,0)Gss{k,0) . (9) 

Using our approximation © for F^^{k,0), approximation (|8a|) for Gnn{k,0) 
and Gssik,0) relations (|lc)) and ()4d|) . one gets after 0-integration: 

Esnik) = qs0oEnik) I Afc , for pn~> p^ , (10a) 

where is given below by ()llb|) . Similarly, in the opposite limiting case: 



Esnik) = qn0oEs{k) / Ak , for p„ » p^ 



(10b) 
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For an arbitrary relation between pn and (including ps ~ pa) we suggest 
the interpolation formula 

EU^) = 0o[qnE,(k) + qsEnik)] / Ak , (11a) 

Afc = ly k"^ + 7nik) + -fs{k) + (qs + qn)0o , (lib) 

which obeys all needed limiting cases. Moreover, in the limiting case 0oqi ^ 
ji it turns into a simpler form 

Esnik) = [psE,{k) + PnEn{k)]/p . (12a) 

This equation has a physically motivated solution 

Esn{k) = Esik) = Enik) , (12b) 

that gives {un{k,t)[un{k,t) — Us{k,t)]) = and thus requires 

Un{r,t) = Us{r,t) , (12c) 

i.e., a fully coherent motion of the superfluid and the normal fluid velocities. 
It means that our interpolation (lllj) is physically justified and actually works 
better than one would expect, having in mind our rather crude approxima- 
tions. 



IC. The energy balance equations 

The energy balance equations can be derived in a standard manner 
from the equations of motion, Eqs. - see e.g. Ref. 11, The only new 
factor in this derivation is the cross- velocity correlations, for which we adopt 
Eqs. ([TT|) . This gives: 

^ + ^ + ..^ (13a) 

= qn0o{^[qsEn{k) + qnEs{k)] - ^n(A:)} 
dEM.dsM 

= qs0o{^[qsEn{k) + qnEsik)] - E,{k)^ , 

with Afc, given by Eq. ()llb|) . The simplest way to model ei{k), suggested 
in Ref. 1161 is to relate Ei{k) and ei{k) in Eq. H13|) in the spirit of the K41 
dimensional reasoning: 

E,{k) = C£,{kf/^k-^l^ . (14) 
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Here C ~ 1 is the Kolmogorov dimensionless constant. In the absence of 
dissipation, Eq. ((Tl)) immediately produces the stationary solution = £ 
with constant energy flux e in the inertial interval of scales. Then Eq. (|14|1 
turns into the Kolmogorov-Obukhov 5/3-law for E'^: 



2. TURBULENT ENERGY SPECTRA IN THE MINIMAL 

MODEL 

2A. Small normal density 

Let us first consider the case 

Ps > Pn , and thus: Qs <^ Qn ■ (16) 

Clearly, in this case the massive superfluid component does not feel the tiny 
superfluid one and thus in the zeroth-order approximation [with respect of 

(Pn/Ps) « 1] 

E,{k) = C7e2/3fc-5/3 ^ j^4i spectrum . (17) 

Indeed, in the limit ^ one has zero in the RHS of Eq. ()13b|) . Then in 
the stationary case dss/dk = 0, i.e.. Eg becomes /^-independent, and from 
Eq. (|14|) one immediately gets Eq. (fTTj) . 

Much more interesting question is about a spectrum of the normal com- 
ponent of the small density, which is essentially affected by the massive su- 
perfluid component. In the limit (jlOj) and for the stationary case Eqs. (|13a() 
and ()llb|) for E^{k) take the form 

<9en(fc) , .2 TP ti,\ ^ l0oqnEs{k) 1 

— uk^ Ery{k) = qj,0Q Er,{k) , (18a) 

ok L Afc J 

Afc = fcV 7n(A;) + is{k) + qn0o , (18b) 
that will be analyzed for the two limiting cases. 

a. Small k =^ full coupling For small k, 7s (being oc k'^^^) is small with 
respect of the (fc-independent) 0oqn- If so, — > Qn0o and 

de^{k)/dk = qn[En{k) - E,{k)] , (19a) 

with the K41 solution 

E,{k) = E^{k)=Ce^/'^k''^/\ = = (19b) 
having full coupling of the velocities, Us = u^. 
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b. Large k ^ decoupling with K41 regime K41 viscous micro-scale, 
kri = defined in the standard manner: 

z.fe2 = 7„(A:,), ^A:,~4/>3/4_ (20) 

In absence of superfluid component, the spectrum of developed turbulence 
of the normal fluid would decay exponentially for k ^ k^i- We demonstrate 
that in the two-fluid system this does not occur due to an additional energy 
flux from the superfluid component to the normal one. Let k > krf. Then 
one simplifies Eqs. (|18|) to an algebraic form with solution 

One finds here a new scale: 

T^kl = gn0o ■ (21b) 

In the case of essential mutual friction, when k^ ^ A;^, (i.e. ^ ^^^)) 
in the subinterval 

k^:$> k:$> k^, (22a) 

a spectrum of normal component (in the former viscous interval!) deviates 
from the K41 just slightly: 

u'^k^ 

E,{k)-E^{k) « ^^^(A;)— (22b) 

It means that in the subinterval (|22aj) the mutual friction is still strong 
enough to keep the normal velocity almost coupled to the superfiuid one: 



k^ 



\Us - -Unl ^ Us ^ n 2 ^ '"s • (22c) 

2gn0o 

However, for 

k>K (23a) 

the strong coupling (|22c() disappears, as it follows from Eq. (|21a() . In this 
region, according to Eq. (|21a() : 

E^^E^^E^^^o,k~^-^l\ (23b) 

i.e. Ea Eg. However, in contrast with standard K41 spectrum, one has 
here a power-law rather than an exponential decay. 
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Fig. 1. The asymptotic forms of the coupled 3D energy spectra (schematic 
log-log plots of energy spectral density versus wave number) in the normal 
fluid (solid black lines) and in the superfluid (solid grey lines) resulting from 
the continuous two-fluid model. The dotted lines indicate the form of original 
conventional uncoupled K41 spectra. The top (bottom) figures correspond 
to limiting cases of large (low) mutual friction, for large (left) and small 
(right) normal fluid density. For further details, see text. 
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2B. Small superfluid density 

Consider now the opposite limiting case to Eq. (fT6|l 

Ps pn, and thus: gs > (Zn • (24) 

Clearly, then (in the zeroth-order approximation) in the inertial interval of 
scales, k < kf], the normal component obeys the standard K41 spectrum 

E^{k) = C7e2/3^-5/3 ^ (25) 

in the inertial interval of scales, k < krj. The spectrum of the superfluid 
component is considered below in various limiting cases. 

a. Small mutual friction Consider first the case when, in addition to 
inequalities (jTO)) . one has 

k^ <.kr^, (26a) 
where the new characteristic scale k^ is defined by 

qs0o = 7s (^t) • (26b) 

For k < k-^ the mutual friction dominate over the nonlinear interaction. 
We already know (see Eq. H12c() ) that in this case = Ug, i.e., there is a 
complete coupling of normal and superfluid velocities, both spectra obeying 
the 5/3-law (IToEl) . 

For k > k'^ the nonlinear interaction dominates over the mutual friction: 

(Us-Mn) « (|«sP) • (27) 

Nevertheless, both spectra obey 5/3-law with Eg — £n- Normal spectrum 
ends at ~ Superfluid component does not feel this cutoff and continues 
further until the limit of the classical description. Beyond this limit, the 
energy cascade is taken over by reconnections and Kelvin waves (see below). 

b. Large mutual friction, In some sense, richer physics corresponds to 
the case with larger mutual friction, when 

A;t > A:^ . (28) 

Then one has a coupled turbulent motion of the normal and superfluid com- 
ponents (with u^{r,t) = Us{r,t)) until the viscous cutoff pO]) . 
For 

kr,<k<k^ (29) 
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one has a case of strong mutual friction with the normal component at rest. 
This case has been considered in details in our Ref. (and agrees with 
numerical result of VinerP^, giving the —3-spectrum 

E,{k) ql0lk-\ (30) 

that originates from the balance of the nonlinear flux and the friction terms 
in Eq. pif) below. 

Indeed, in this case one puts in Eq. ()13b|) Ea = 0, ~ 0o'7s, and the 
RHS of Eq. H13b|) can be approximated as 

RHS ^ 0oqsEj- - l) ^ -0oqsEs . 



In that case, instead of Eq. H13b|) . one arrives at 

qs0oEs{k) = ^ = C-y^ ^ [k^/'E-i/Hk)] , (31) 

which has the solution H30|) . 

For k > k-^ the mutual friction is NOT important and the superfluid 
component recovers the 5/3-law, but with smaller energy flux. 



3. ROLE OF KELVIN WAVES 

The description we presented above ignores the fact that turbulence of 
the superfluid component consist of discrete tangled vortices with quantised 
circulation k. This description is valid for the scales greater that the mean 
distance £ separating the quantised vortices. For example, in '^He above 
1.5 K the normal fluid Kolmogorov scale appears to be of the same order 
as i and, therefore, in this case our model is not applicable for predicting 
the superfluid spectrum below the Kolmogorov scale^^. Another example 
is turbulence at T < IK when the normal component is extremely weak 
and behaves like a Knudsen gas rather than a fluid. In this case, most of 
energy reaches £ along the cascade without dissipation and one should ask 
what happens to this energy below this scale. At scale £, an essential role 
in turbulence evolution play vortex reconnections during which part of the 
energy is lost to phonon emission, and the rest of the energy is transferred 
to Kelvin waves, see e.g. paper of W.F. Vinen in this issue. The recon- 
nections produce sharp cusps which quickly transform into a superposition 
of Kelvin waves whose nonlinear interaction leads to further turbulent cas- 
cades through scales. Note that these sharp cusps correspond to a broad 
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distribution in wavenumber space and, therefore, both the direct energy cas- 
cade and the inverse cascade of waveaction can be important in subsequent 
evolutior^. To describe the statistical nonlinear Kelvin waves one can use 
the weak turbulence approach, which results in a six-wave kinetic equation 
for the energy spectrun^. A differential equation model for Kelvin wave 
turbulence (kelvulence) which preserved the essential scalings and solutions 
of the original integral kinetic equation was derived in Ref 1171 

where n = EgL/to is the waveaction spectrum^, L is the vortex length per 
unit volume - the vortex line density, k is the circulation quantum, C is a 
dimensionless constant and uj = uj{k) = KA;^/(47r) is the Kelvin wave fre- 
quency^ . This equation preserves the energy (per unit length of the vortex) 

E = f oj'^/^nduj (33) 

2v n J 

and the waveaction (per unit length of the vortex) 

= I uj'^l'^ndL) . (34) 

2^ J ^ ' 

Equation (|32)) has both the direct cascade solution n ~ k~^'^^^ and the 
inverse cascade solution n ~ k~^. It also has the family of thermodynamic 
Ray leigh- Jeans solutions, 

T 

n = . (35) 

u; + M 

where T and ^ are constants having a meaning of temperature and the 
chemical potential, respectively. 

The direct cascade of energy in kelvulence is eventually dissipated at 
high wavenumbers either via phonon radiation or via friction with the normal 
component. The phonon radiation is the dominant dissipation mechanism 
near absolute zero whereas the mutual friction becomes more important at 
higher temperatures, e.g., at T > 0.4 K in ^He, as estimated in Ref^ Using a 
dimensional argument and an assumption that the radiation is quadrupolar, 
we have the following expression for the sound dissipation terir^, 

in)rad = TTTT' (^^^ 

^Factor L appears when one calculates energy per unit volume in terms of the energy 
per unit vortex length 

■^Here, we ignore logarithmic factors 
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where Cs is the speed of sound. This kind of sound radiation terminates the 
energy cascade at a finite frequency Wrad ~ {e^cf / k^^)^^^^ , where e is the 
total energy injection rate per unit vortex lengthP^. 

Now let us introduce the effect of mutual friction. It was argued in 
Ref n that the characteristic time of the mutual friction dissipation at low 
temperatures, when the normal fluid is at rest, is ~ l/{aKk'^) where a is 
the (temperature dependenlP^H) mutual friction coefficient. Thus, we can 
postulate the following dissipation term in our differential model, 

{n)rai = —ank'^n. (37) 



Interestingly, this term looks like a viscous dissipation with an effective 
viscosity coefficient an. Its temperature dependence qualitatively agrees 
(within a numerical factor of order unity) with measurements of the effec- 
tive kinematic viscosity in ^He extracted from the towed grid experiments in 
the range between about 1.2 K to 1.7 KP^. For higher temperatures, coeffi- 
cient a has to be replaced with a more complicated function of both friction 
parameters a and a' For a finite counterflow velocity Vns, one should 
replace ^ witlJ^U 



(n) 



mf 



2lo g(l/fca) 
47r 



n . 



(38) 



which describes Glaberson instability of Kelvin waves when their phase ve- 
locity is less than V^s- Finally, we are leaving for future consideration an 
interesting case of turbulent Glaberson amplification when V^s is random 
due to turbulence in the normal fluid. Note that the friction dissipation in 
()38() differs from 1)3 7|) by factor log(l/A;a)/47r which, following Ref^ we will 
assume to be close to unity. 

To examine the effect of mutual friction on the energy cascade, let us, 
following Ref. 1171 introduce a reduced second order model which ignores the 
waveaction conservation and the inverse cascade. With the friction term 
()37() included, we have 



Ci d 



n 



..10 



duj 



(39) 



where Ci > is a dimensionless constant. The general stationary solution 
of equation (|39|) is 



n 



2e 

c7 



4/5 



,42/5 



47raK^° 



4/5 



1/4 



1 

17/10 



(40) 
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We see that at low frequencies this expression coincides with the non-dissipative 
energy cascade spectrum, n ~ uj~^'^^^^, and we also see a sharp cut-off at 

= (4vra)V4.2 • (41) 

Comparing ujrad and one can see that the frictional dissipation becomes 
more important than the phonon radiation if 

e < (4^a)i3/8^2^ ^ (42) 

This expression differs from expression (87) of Ref. 1 (assuming their relation 
(74) between e and I, they obtained crossover between the two regimes at a ~ 
7 X 10^*^ which corresponds to the temperature about 0.4 K). We attribute 
this difference to the fact that Ref. Q assumed the sound radiation to be 
dipolar rather than quadrupolar as in Ref. [T71 

Here we present our estimate of the He H temperature when condition 
(|42j) indicates a crossover between phonon radiation and mutual friction 
dissipation: injecting a power at the level of 1 W into 1 liter of liquid of 
density about 145 kg/m^ (i.e., energy decay rate e 7 m^/s^) results to 
a typical vortex line density 10^'' — 10^^ m~^, giving e ~ 10^^'^ m^/s^. 
The condition (|42|) thus requires Ana of order 10~^ (i-e., about an order of 
magnitude higher than the estimate of Ref. |T]) , which roughly corresponds^ 
to temperature of about 0.5 K . 



4. CONCLUSIONS 

In this paper, we presented a minimal model for turbulence of the cou- 
pled superfluid and normal components in superfluid helium. The model 
comprises a system of nonlinear partial differential equations for the energy 
spectra and its origin in the case of classical fluids can be traced back to Ko- 
vasznay 1947 papei^. The basic idea of such models is that the nonlinear 
terms, being of the simplest possible form, should preserve the original turbu- 
lence scalings and, in particular, predict correctly the Kolmogorov cascade. 
Having the Kolmogorov scalings built into the model and adding additional 
physical interactions, such as mutual friction and viscosity, one then obtains 
new nontrivial physical regimes charactrerised by non-Kolmogorov spectra. 
For superfluids, the first application of such model was done in Ref. 1111 in 
the limiting case of the normal fluid at rest, and the new —3 spectrum was 
predicted. 

In the present paper, we generalised this model to the case where both 
the normal and the superfluid components may be turbulent. The crucial 
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theoretical step here is our estimation of the cross-correlation function be- 
tween normal and superfluid velocities which determines their joint dynam- 
ics. This cross-correlations leads to appearance of a rich variety of interesting 
new regimes. In particular, we found for large superfluid to normal density 
ratio a regime in which the normal component is "dragged" by the superfluid 
component via the mutual friction with resulting Kolmogorov spectrum ex- 
tending far below the Kolmogorov dissipation scale. Our model also confirms 
the picture suggested in Ref . jl2j of a "knee" spectrum where the dissipative 
—3 scaling at medium wavenumbers exists in between of —5/3 Kolmogorov 
ranges at low and large fc's. Our theory bridges classical turbulence with 
quantum turbulence and in a quantitative manner, it points out similarities 
and differences between the two, and we expect it to be useful and efficient 
for numerical simulations of more complicated experimental cases. 

We also discuss the case when our continuous two-fluid description 
breaks down at scales below the mean distance between the quantised vortex 
filaments. At these scales, the turbulent cascades are believed to be carried 
through by random nonlinear ly interacting Kelvin waves (kelvulence) . A 
differential model for kelvulence, including the phonon radiation effect, was 
proposed in Ref. '171 In the present paper we extend this model to include the 
mutual friction effect and we obtain an analytical solution where the energy 
cascade is arrested at a finite wavenumber. Comparing this friction cut-off 
with the previously obtained radiation cutofP^ we obtained an estimate for 
the crossover between the radiation and the friction dissipation mechanisms. 

We leave the exact details of applicability to He II and to '^He-B for 
future work, as well as an extension of our approach to (anisotropic) coun- 
terflow turbulence. Let us point out that the likely candidate of experimental 
verification of our theory might be ^He-'^He superfluid mixtures where, due 
to presence of ^He quasiparticles, mutual friction at low temperature is still 
expected to be significantly higher than in pure He II. 
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